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The  stability  of  cytosine  in  aqueous  solution  was  studied  in  the  laboratory,  simulating  prebiotic 
conditions  and  using  gamma  radiation  as  an  energy  source,  to  describe  cytosine  behavior  under 
radiation.  For  a  better  understanding  of  the  radiation-induced  processes,  we  proposed  a  mathematical 
model  that  considers  chemical  reactions  as  nonlinear  ordinary  differential  equations.  The  radiolysis 
can  be  computationally  simulated  by  an  agent-based  model,  wherein  each  chemical  species  involved  is 
considered  to  be  an  agent  that  can  interact  with  other  species  with  known  reaction  rates.  The  radiation  is 
contemplated  as  a  factor  that  promotes  product  formation/destruction,  and  the  temperature  determines 
the  diffusion  speed  of  the  agents.  With  this  model,  we  reproduce  the  changes  in  cytosine  concentration 
obtained  in  the  laboratory  under  different  irradiation  conditions. 
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1.  Introduction 

Radiation-induced  oxidation  reactions  are  a  crucial  process 
in  understanding  the  formation  of  biologically  relevant 
molecules  in  planets  and  icy  bodies  such  as  comets  [1- 
4].  Primitive  Earth  and  astrophysical  conditions  for  these 
processes  are  difficult  to  simulate  in  the  laboratory.  There  are 
few  experimental  setups  studying  the  stability  and  formation 
of  some  compounds  of  biological  importance  under  a  high 
radiation  field  [5,  6].  In  these  experiments,  changes  induced 
by  ionizing  radiation  can  be  quantified  by  dose  measurements 
that  indicate  the  amount  of  energy  deposited  on  the  samples 
by  the  gamma  radiation.  The  first  experiment  of  prebiotic 
synthesis  was  proposed  by  Miller  [7],  where  several  amino 
acids  were  synthesized  from  different  gases  in  a  primitive 
reducing  atmosphere  (H20,  CH4,  NH3,  H2)  and  electric 
discharges.  Prebiotic  importance  for  molecules  specially 
RNA  was  established  by  Gilbert  [8].  One  component  of 
RNA  and  DNA  is  cytosine,  a  pyrimidine  base,  but  there  are 
few  experiments  about  the  stability  of  cytosine  (C4H5N30) 
under  radiation,  the  molecule  of  interest  in  this  work. 

The  irradiation  of  cytosine  in  aqueous  solution  involves 
the  interaction  of  cytosine  molecules  with  the  different  species 
formed  by  the  radiation-induced  decomposition  of  water  (H, 
OH,  ea“,  H2,  and  H202).  To  describe  the  products  generated 


by  the  interaction  of  the  different  species  under  radiation  in 
an  aqueous  medium,  we  propose  a  mathematical  model  that 
describes  the  mass  balance  of  all  species  involved,  considering 
chemical  reactions  as  nonlinear  ordinary  differential  equations 
(NODE)  [6,  9].  This  model  is  complicated  due  to  the 
significant  number  of  reactions  involved,  the  coupling  between 
equations,  and  by  the  fact  that  all  the  NODEs  need  to  be 
solved  simultaneously.  Moreover,  the  non-linear  character  of 
the  equations  makes  the  model  strongly  dependent  on  initial 
conditions.  To  circumvent  these  issues,  some  authors  have  used 
Monte-Carlo  simulations  [10].  We  have  proposed  an  agent- 
based  model  [11,  12]  to  simulate  the  chemical  evolution  of 
oxidation  reactions  of  ferrous  ions  under  radiation.  The  model 
is  a  modified  version  of  the  prey-predator  model  [13,  14]  in 
which  the  mass-balance  equation  includes  sink  terms  (all  the 
reactions  leading  to  destruction  that  can  be  considered  prey) 
and  source  terms  (all  the  reaction  rates  leading  to  production 
that  can  be  interpreted  as  predators).  In  this  model,  each 
chemical  species  involved  is  considered  as  an  agent  that  can 
interact  with  other  species  with  known  reaction  rates,  radiation 
is  taken  as  a  factor  that  promotes  a  products  formation/ 
destruction,  and  the  temperature  determines  the  diffusion 
speed  of  the  agents.  Here,  we  modify  the  model  to  reproduce 
the  radiation-induced  reaction  of  cytosine  in  aqueous  solution. 
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2.  Experimental  Setup 

The  experimental  setup  followed  the  setups  of  [14, 
15].  Aqueous  solutions  of  5  X  10“4  M  of  cytosine  were 
irradiated  with  a  Cobalt  60  source  at  an  irradiation  dose 
ranging  from  1  to  5  kGy.  The  samples  were  in  Pyrex  cells  at 
room  temperature  (295  K),  oxygen-free.  After  irradiation, 
the  samples  were  analyzed  by  UV-Vis  spectroscopy  at 
275  nm  (the  peak  wavelength  associated  with  cytosine 
[16])  in  a  Varian  Caryl  00  spectrophotometer  (California 
USA).  Cytosine  concentrations  were  determined  and 
compared  with  a  modified  version  of  the  agent-based  model 
programmed  with  Python  in  [11,  12]. 

3.  Model  for  Radiation  Induced  Reaction  of 
Cytosine 

3. 1  Chemical  Reaction  Model 

Chemical  reactions  involved  on  cytosine  production  involve 
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The  chemical  reactions  involved  in  our  cytosine  experiment 
starts  with  a  reaction  that  can  occur  with  a  probability  of 


87%: 

[1]  +  [2]A[3]  (1*) 

or  a  reaction  that  occurs  with  a  probability  of  10%: 
[l]  +  [2]-[4]  (1*) 

Then  the  following  reactions  occur: 

[3]  +  [2]4[5]  (lr) 

[4]  +  [2]  ->  [5]  (1  d) 

[l]  +  [6]2[7]  (1.) 

[8]  +  [9]2[l0]  (1)5 


Now  the  next  reaction  can  occur  with  a  probability  of  82%: 
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[10]  +  [2]4[11]  (lj) 

or  the  reaction  with  a  probability  of  18%: 

[10]  +  [2]4[12]  (\h) 

Afterwards  the  next  reactions  occur: 

[11]  +  [2]4[13]  (If) 

[12]  +  [2]4[13]  (1 j) 

with  the  six  known  reaction  rates: 


number  6.022  X  1023  mol,  k  is  the  Boltzmann  constant 
1.38  X  10-23  m2kgs_2K_1,  and  G.  is  the  radiochemical 
constant  that  gives  the  number  of  species  (/)  produced  when 
1 00  eV  are  absorbed  by  the  system.  For  the  case  of  cytosine, 
G  is  4.42  molecules,  then: 

/,(/„  T)  =  [7.01-  2.04/,  ]xlO-4M,.,  (4b) 

Chemical  reactions  described  by  the  system  of  Eqs.  (1)  can 
be  rewritten  using  Eqs.  (3)  and  (4)  as  the  following  NODE 
system: 

^l  =  /,-r,[l][2]-,,[l][6].  (5„) 


j •  =  6.3  x10s, 

(2d) 

r2  =  3.1x10s, 

(2b) 

r3  =  13x10s, 

(2c) 

r5  =  0.47x10s, 

(2d) 

r6=  5.7  x10s, 

(2c) 

r7  =3.0xl09. 

(2 f 

-r7|ll][2]-fl[12][2], 


^-/,+M  1][2]-<1[3][2], 

Of 

^-/.+-l[l][2]-r,[4l[2]. 

0  d) 

^  =  /,+<!  M2]  +  ^[4][2], 

Oe) 

3.2  Mathematical  Model 

Kinetics  of  the  radiation  induced  chemical  reactions  in  a 
thermal  bath  (Eqs.  1)  involves  the  mass  balance  for  each 
of  the  13  species.  This  balance,  assuming  that  is  due  with 
the  interaction  between  pairs,  can  be  described  by  coupled 
nonlinear  differential  equations  considering  external  source 
(radiation  and  bath  temperature),  production  (source  terms) 
and  destruction  (sink  terms),  as  proposed  on  a  previous 
work  [6]: 


m 

dt 

d[ 7] 

dt 

m 

dt 


=  f6-rAl¥\+r2[^¥\’ 
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=  fs~r5[  8]  [9], 
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where  N  denotes  the  molar  concentration  of  the  i-species 
at  time  t,  r(i{  denotes  the  reaction  rate  constant  for  the  two 

j>k  ,  .  . 


external  source  given  by: 

/■(/,, r)=6-2Xl°"  G,e-  M 
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molecular  mass  of  water  (18.02  g/mol),  NA  is  the  Avogadro 
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Figure  1.  Cytosine  concentration  as  function  of  the  radiation  doses. 


3.3  Agent-Based  Model 

The  coupled  system  of  Eqs.  (5)  can  be  solved  using  an  agent- 
based  model  like  the  prey-predator  one  that  we  developed 
in  previous  works  [11,  12]  where  we  generated  a  Python 
code  to  evaluate  the  molar  concentrations  of  each  species. 
In  Fig.  1  is  plotted  the  cytosine  concentration  obtained  by 
the  code  comparing  the  values  with  the  experimental  ones 
for  different  radiation  doses.  As  we  see  our  numerical  result 
closely  resembles  the  experimental  one,  with  a  difference  of 
less  than  10%.  Red  stars  are  concentrations  from  the  agent- 
based  model  and  black  dots  from  experimental  data. 

4.  Conclusion 

Cytosine  reactions  induced  by  gamma  radiation  at 
room  temperature  were  studied  and  compared  with  the 
concentrations  determined  by  an  agent-based  model.  The 
model  evaluates  the  concentration  of  products  generated 
by  the  interaction  of  different  reactive  free-radicals  under 
radiation.  It  involves  the  mass  balance  of  16  species, 
considering  each  species  as  an  agent  that  can  interact  with 
other  species  with  known  reaction  rates.  Interactions  lead 
to  destruction  (prey)  and  production  (predator)  terms, 
with  radiation  considered  as  a  factor  that  affects  product 
formation.  This  simple  and  robust  model  agrees  with  the 
experimental  results.  Both  approaches  showed  that  cytosine 


decomposed  rapidly  in  a  high  radiation  field  environment, 
and  that  for  the  survival  of  this  molecule  and  its  further 
participation  in  the  formation  of  more  complex  molecules, 
it  is  necessary  to  have  a  protection  mechanism,  such  as  the 
adsorption  in  clay  minerals. 
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